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Studies aimed at understanding the global properties of the hyperpolarizabilities have focused on 
identifying universal properties when the hyperpolarizabilities are at the fundamental limit. These 
studies have taken two complimentary approaches: (1) Monte Carlo techniques that statistically 
probe the full parameter space of the Schrodinger Equation using the sum rules as a constraint; 
and, (2) numerical optimization studies of the first and second hyperpolarizability where models of 
the scalar and vector potentials are parameterized and the optimized parameters determined, from 
which universal properties are investigated. Here, we employ an energy spectrum constraint on the 
Monte Carlo method to bridge the divide between these two approaches. The results suggest an 
explanation for the origin of the factor of 20-30 gap between the best molecules and the fundamental 
limits and establishes the basis for the three-level ansatz. 
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I. INTRODUCTION 

The fundamental limit of the off-resonant electronic 
hyperpolarizability and second hyperpolarizability are 
calculated using the constraints on the energies and ma- 
trix elements of position operator (which we loosely call 
them transition moments throughout this text) imposed 
by the Schrodinger Equation. [l-Q It is most convenient 
to quantify these constraints in the form of the Thomas- 
Kuhn sum rules, and then using the assumption that 
when first and second hyperpolarizabilities, /3 and 7 re- 
spectively, are optimized only three states (including the 
ground state) contribute - called the three-level ansatz. 
The upper limit of /3 is given by, 
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where — e is charge of the electron, h is Planck's constant, 
N is the number of electrons of the quantum system and 
Eio is the energy difference between the first excited state 
and ground state. The second hyperpolarizability follows 
along the same lines as for the hyperpolarizability, [l|, 0, 
13, m and yields. 
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A comparison of the largest experimentally measured 
hyperpolarizability with the fundamental limit reveals a 
large gap between the two. Measurements have never 
crossed the limit and are typically well below it. 0-0] Un- 
til about 2007, the hyperpolarizabilities of all molecules 
ever measured mysteriously fell a factor of about 30 
below the limit. As such, the limits were seen to 
provide an absolute metric of the nonlinear-optical re- 
sponse. This idea was used, for example, by Slepkov and 



coworkers to understand the surprising nonlinear opti- 
cal properties of polyene oligomers; [1] and by May and 
coworkers to show the promise of small molecules. [1, 
Chen applied the limits to study the nonlinear re- 
sponse of nano-engineered polymers. [ll| More recently, 
a new class of twisted molecules were reported by Kang 
and coworkers that appear to have exceptionally large 
hyperpolarizabilities. |12l4l5l | 

A more important consequence of the limits is the fact 
that they provide a method for determining how the 
nonlinear-optical response scales with the quantum size 
of the system. It is straightforward to show that if the 
scalar potentials and energies are re-scaled by a factor 
fi^^ and the vector potential and positions are all rescaled 
by a factor fi, the shape of the wavefunction remains the 
same aside from it being compressed by the scaling factor 
/i.[T6| This transformation leaves the intrinsic hyperpo- 
larizabilities, defined by. 
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invariant. As such, the intrinsic hyperpolarizabilities can 
be used to compare molecules of very different shapes 
and sizes. 

The concept of scale invariance suggests that opti- 
mization of the hyperpolarizability is a two-step process. 
First, a successful paradigm is identified that optimizes 
the intrinsic hyperpolarizability. Subsequently, a quan- 
tum system within that paradigm can be made larger us- 
ing simple scaling. As an example, Roberts and cowork- 
ers have shown that triphenylamine-cored alkynylruthe- 
nium dendrimers have a nonlinear-optical response - as 
characterized by two-photon absorption (TPA) cross sec- 
tion per electron - that is about an order of magnitude 
larger than the non-triplienylamine-cored versions. [17 
However, when the results are properly scaled according 
to the intrinsic two-photon cross-section, [isi [l9| all of the 
triphenylamine-cored dendrimers had the same intrinsic 
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value that was found to be two orders of magnitude larger 
than the non-triphenylamine-cored versions, poj Thus, 
the triphenylamine-cored dendrimer is the new molecular 
paradigm with a larger intrinsic TPA cross-section than 
other systems; and, dendrimers with triphenylamine- 
cores can then be made larger by adding additional den- 
drimer generations (i.e. additional branches) to increase 
the absolute TPA cross-section. 

Scale invariance has defined two approaches for build- 
ing an understanding of the nature of a quantum sys- 
tem when the hyperpolarizability is near the limit, 
both of which seek to optimize the intrinsic quantities. 
The first are studies that use a parametrization of the 
scalar and vector potentials to determine the nature 
of the potentials that lead to an optimized nonlinear 
response. |2l| This approach has been used to suggest 
a new paradigm for making molecules [12, that has 
lead to the synthesis and characterization of a record 
intrinsic hyperpolarizability. [l^, Such studies also 
show that a quantum system near the fundamental limit 
shares certain universal properties. [2^ Wang and cowork- 
ers have used related optimization techniques to theo- 
retically build larger molecules from smaller molecular 
building blocks as the basic units. [27| 

A second approach is using Monte Carlo simulations 
that statistically probe the full space of allowed solutions 
to the Schrodingcr Equation by using the Thomas-Kuhn 
sum rules as a constraint. (28l 129} In contrast to the po- 
tential function approach, which yields an optimized in- 
trinsic hyperpolarizability of 0.78, the Monte Carlo ap- 
proach yields a distribution of hyperpolarizabilities that 
approaches arbitrarily close to the fundamental limit. 
The Monte Carlo approach is the most general, but, gives 
no information of how to design a quantum system with 
a large hyperpolarizability. Our present work seeks to 
bridge the divide between the Monte Carlo approach and 
real systems by classifying materials according to their 
energy spectrum. 

In contrast to studying specific system to gain an 
understanding of the origins of the nonlinear optical 
response. [30l432l | our work seeks to identify broad prin- 
ciples with the aim of applying these results to finding 
new paradigms of the nonlinear optical response. [13 



rules, which are generalizations of Equation |4l [SSHSSj 
are not required for most molecular systems. Sum rules 
have also been used to study the dispersion of the non- 
linear response; [3^ but, our present interest is in the off- 
resonant regime. 

For simplicity, we use the dimensionless sum rules 



E 



(5) 



where = Eio/Eio (the ground state is labeled by 0) 
and the normalized transition moments are given by 
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Equation!?] defines the upper bound of the transition mo- 
ment from the ground state, xqi, by use of Equation |4] in 
light of the fact that Eiq < Eio for j > 1. 

The numerical procedure we apply here is similar to the 
previous studies of the hyperpolarizability, [2^ and the 
second hyperpolarizability, [29| with the exception that 
the energy levels are not chosen randomly. Instead, we 
set the energies to be of a predefined functional form. 
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For a particular function /(s), the Monte Carlo simula- 
tions span a restricted set of solutions of the Schrodingcr 
Equation that gives a statistical sampling of the solutions 
that have a particular energy spacing. This approach of 
fixing the energy-level spacing is motivated by the obser- 
vation that detuning in real systems from optimal energy 
spacing is responsible for the factor of 30 gap between the 
best molecules and the fundamental limit, [y, 0, 113 

The energy values in our model must be chosen for 
states j > i such that Ej > Ei. Thus, for energy classes 



such as Ei 



we shift all energies so that this 



constraint holds. As an example. 



II. APPROACH 

The Thomas-Kuhn sum rules, which are a direct con- 
sequence of the Schrodingcr Equation, are given by. 
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where 5ij is the Kronecker delta, Em is the energy of 
state m and x„,„ is the (n, m) position matrix element. 
Equation|3]holds for multi-electron Hamiltonians with N 
electrons with spin, externally app lied electromagnetic 
fields, and electron correlations. [I3 The relativistic sum 
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Thus, the ground state is labeled by j = and has zero 
energy. Since = Eiq/ Eiq, then cq = and ei = 1. 

The transition moments are assigned, as in our pre- 
vious work [2^, [2^ , so that the energies and transition 
moments are consistent with sum rules. Since we use 
the dipole-free expression for /3 and 7 [s^ we need 
only use diagonal sum rules (to = p) in Equation [S] to 
get all of the required transition moments. Starting with 
(to,p) = (0,0), we get 

eio|ao|' + e2o|6o|' + e3o|e3o|' + ... + e„o|?«o|' = 1- (10) 



eio = ei — eo = 1, so ^oi is randomly assigned from the 
interval —1 < ^oi < 1- Subsequently, ^02 is obtained 
from 
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do < (1 - eiCi'o) /e20. 
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A random number —1 < r < 1 is used to get ^20 from 
Equation [121 
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E = exp(j) 
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For a system with s-states, the procedure is repeated for 
all other transition moments, ^io, except for which 
is directly determined from the last remaining diagonal 
sum rule. 

The next set of transition moments are calculated us- 
ing the higher-order sum rules in sequence with p = 
1, 2, 3, . . . in Equation [5l Since is assumed to be real 
for all states i and j. 
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The energy and transition moment values are then used 
to calculate the intrinsic first and second hyperpolariz- 
abilities. This procedure is iterated repeatedly to gener- 
ate a distribution of first and second hyperpolarizability 
values. The process is then repeated for a variety of en- 
ergy level functions. 

Using the dipole-free sum over state expression [sl] , (3 
in the off-resonant regime is given by 
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where a prime indicates the ground state is excluded 
from the summation. The second hyperpolarizability, 7 
is given by 
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FIG. 1. Distribution of Pint for energy function Ej oc exp(j). 



III. RESULTS AND DISCUSSION 



We apply Monte Carlo simulations to eight classes of 
energies, -j"^, j^^^, j, f, and The 

nature of their distributions can be used to determine 
what properties are relevant for optimizing the nonlinear 
response. 



A. The Three-Level Ansatz 

The three-level ansatz states that when the nonlin- 
ear response of a quantum system is at the fundamen- 
tal limit, only three states contribute. fl6l| It can be 
rigourously shown that the linear response is optimized 
for a two-level system. However, for a nonlinear response, 
an analysis of the two-level model shows that it violates 
the sum rules. The three-level ansatz was used in calcu- 
lating the fundamental limits based on the argument that 
an optimized system concentrates the oscillator strength 
in the least number of states. Since the three-level model 
is the lowest-order approximation that obeys the sum 
rules, it was used without proof. [J 

While the three-level ansatz has not been proven an- 
alytically, it appears to always hold in a large number 
of simulations. We find the same results here. Figure 
[T] shows the distribution of /? for an exponential energy 
function Ej = e-* . Shown are runs for 3 to 80 states where 
the distributions for 20 and 80 states overlap. Thus, 
adding more states will not change the character of the 
distribution. The three-level model shows a distribution 
that peaks at the limits while all of the other distributions 
peak at Pint = in a cycloid-like function, a result that 
we observed in previous studies. Interestingly, all these 
distributions tail off dramatically at /3int = 1- These 
results are clearly consistent with the three-level ansatz. 

Figure [5] shows the distribution of "fint , also for an ex- 
ponential energy function Ej = . Note that the second 
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FIG. 2. Distribution of 7i„f for energy function Ej oc exp(j). 

hyperpolarizability has a different positive and negative 
limit. As in the case of /?, the distribution of the three- 
level model peaks at the fundamental limit, where the 
second hyperpolarizability is positive. Again, all of the 
other distributions peak at 7i„t = 0. Thus, the three- 
level ansatz also holds for the second hyperpolarizability. 

The shape of the distribution changes dramatically 
near = 1 and jmt = 1 as a function of the num- 
ber of states. For n — 3, the function is strongly upward 
sloping even on the log plot, showing that statistically, 
a three-level system with exponential energy spacing is 
most likely to be found near the fundamental limit. In the 
four-level model, the distribution is flat over most of the 
domain, suggesting that all values of the intrinsic hyper- 
polarizabilities are equally likely, but is strongly sloped 
downward near the limit. As the number of states is in- 
creased, the distribution becomes more strongly peaked 
near zero first and second hyperpolarizabilities. This be- 
havior shows how a many-level system will almost always 
have a smaller intrinsic first and second hyperpolarizabil- 
ities than a three-level system. 



B. Hyperpolarizability 

Analytical calculations of the fundamental limits using 
the three-level ansatz suggest that the ideal system is one 
in which the ratio E10/E20 — > 0, i.e. when the second 
excited state energy is much larger than the first. This 
type of energy-level spacing is not typically found in real 
quantum systems. For a particle in an infinite square 
well, the energy scales as Ej — and for a harmonic 
oscillator it scales as Ej = j. For a single electron and 
point nucleus, on the other hand, the energy scales as 
Ej ~ j^"^ ■ We define an energy class by the functional 
form of the energy level spacing. Thus, the harmonic 
oscillator falls in the energy class Ej ~ j and one electron 
atoms in Ej = ■ 

Calculation of the fundamental limit of the hyperpo- 



larizability leads to [l|, , 

P - 6^3^^ ' ^^2^ ' G{X)f{E) = p^G{X)f{E), (17) 
where ^0 = fiynax, 

f{E)^{l~Ef/^{E^ + ^-E + iy (18) 

and 

G{X) = ^^^^(1-^4), (19) 

with X = xiqIx^q^^ and E = Eio/E2a- At the limit, 
G{X) = 1 and f{E) = 1. 

Each energy class corresponds to a specific value of E 
and therefore a specific value of f{E). In the case of the 
three-level model, each random assignment of transition 
moments simply identifies one value of X. For example, 
for the exponential energy class, E = 0.27 and f{E) ~ 
0.922. The distribution of intrinsic hyperpolarizabilities 
asymptotically approaches /3i„t = 0.922 for those cases 
where G{X) = 1. For the three-level model, the largest 
value of the intrinsic hyperpolarizability for a particular 
energy class is thus /3i„t = P/Pmax = f{E). Thus, f{E) 
can be interpreted as the largest intrinsic value that is 
allowed by the three- level ansatz. 

Figure |3] summarizes the Monte Carlo simulations of 
several energy classes. The dashed vertical lines represent 
Pint = if{E). The three-level model is thus observed to 
behave as predicted - the ceiling in the intrinsic hyper- 
polarizability gets larger as E increases. 

In general, an n-level model will have n — 2 
parameters. [23j Thus, one would expect that with more 
parameters, it may be possible to find combinations that 
yield a larger intrinsic hyperpolarizability than the three- 
level model. Interestingly, the three-level model yields a 
distribution that peaks at Pint = if{E). Thus, if the 
energy spacing is optimized, the distribution peaks at 
Pint ~ ±1- In contrast, when more states are included, 
the distribution peaks at Pint = and falls off sharply 
at Pint = ^f{E)- Thus, while the energy function f{E) 
was specifically defined for the three-level model, it ap- 
pears to define the cutoff for any number of levels for any 
quantum system that obeys the sum rules. 

There are practical implications of these observations. 
First, when a quantum system is designed to concentrate 
oscillator strength in only two excited states, from a sta- 
tistical perspective, it is most likely that the resulting hy- 
perpolarizability will be near the limit. If the oscillator 
strength is shared amongst many states, then the most 
likely hyperpolarizability is zero. Since molecules typ- 
ically have many states with large transition moments, 
it is not surprising that the hyperpolarizabilities fall far 
short of the fundamental limit. Note that beyond 20 
states, all the resulting distributions appear identical. 

When the energy function is suboptimal, the intrinsic 
hyperpolarizability can exceed the value obtained from 
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FIG. 3. Distribution of Pint as a function of the number 
of states for several different energy functions. The vertical 
dashed lines represent f{E) for that energy function. 



TABLE I. Summary of the properties of the hyperpoiarizabil- 
ity for various energy functions when the hyperpoiarizability 
is small or near the limit. 



Function 


Pint 


Levels 


Dominant 
States 


Maximum 

rtint 
Pnm 


exp(j) 


-0.981 


3 


1,3 


Pif = -0.78 


j 


-0.90 


4 


1,8,9 


Pit = -0.56 




0.284 


9 


1,4,5,6,7,8,11,12 Pit = -0.310 



a three-level model, that is, Pint > f{E). This suggests 
that it may be possible to make a quantum system with 
a large intrinsic hyperpoiarizability in which many states 
contribute. It is therefore interesting to investigate the 
proprieties of the outliers in the distribution to determine 
if this may lead to a new paradigm for making molecules 
with a large hyperpoiarizability. 

Table |I] summarizes the properties of an outlier for 
the energy class exp(j), in which the hyperpoiarizabil- 
ity exceeds the three-level model [Pint = —0.981 and 
|/(£;)| = 0.922), the energy class j with |/(£;)| = 0.707 
and class of energy a case were the intrinsic hy- 
perpoiarizability is much less than the for a three-level 




FIG. 4. PtS^ for the energy class j"^ when Pint = 0.284. 

system. We use /3nm, a single term in the double sum in 
Ecj^uation I15[ to represent the fractional contribution to 
the hyperpoiarizability of the two states n and m. The 
total hyperpoiarizability is then given by, 

/3 = 5I/5"™- (20) 

n,m 

Interestingly, states 1 and 3 are responsible for 80% of 
the hyperpoiarizability. Thus, while this outlier was not 
constrained to be a three-level system, three states dom- 
inate the response. In contrast, when the hyperpoiariz- 
ability is well below the limit, many states contribute. 
The second entry in Table U shows that 8 excited states 
contribute to yield (3int = 0.284 for the energy class 
for which \f{E)\ = 0.184. FigureHshows a plot of 
In this case, the dominant contribution is of opposite sign 
and of larger magnitude than of the total nonlinear re- 
sponse. This illustrates how many states may contribute, 
but in such cases, the contributions can cancel. 

Figure [5] shows a density plot of the hyperpoiarizabil- 
ity as a function of xqi/xqI""^) where xj}^"^ is the 
fundamental limit of the transition moment defined in 
Equation [7l For each energy class, the Monte Carlo sim- 
ulations generate all possible values of X in the [—1,1] 
range as well as all the possible transition moments be- 
tween all pairs of states. Thus, the figure represents a 
sampling of the subset of the full Hilbert space for that 
energy class. Associated with each energy class is a sin- 
gle value of E. The plots in order of etc. 
(left to right, top to bottom) are in order of decreasing 
E. 

We can interpret this behavior by combining Equations 
[TTlfTOl which come from the three- level ansatz, to get 

/3»t = G{X)f{E) 

= mi - Ef'^ (e^ + + 1) 

(21) 
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FIG. 5. Density plots from 1 million Monte Carlo runs of 
Pint as a function of X for several energy classes of a 20-state 
model. The dashed lines represent ±f{E). Red, blue and 
green colors correspond to the Pint values for which 3, 4 to 6 
and 7+ states dominate the nonlinear response, respectively. 
The black solid curve represents a plot of ±/(_E)G(X) for E 
fixed by the energy class and as a function of X. 



It is worthwhile to dwell on the origin of Equation [211 It 
is derived from the three-level model of f3 using the sum 
rules to reduce the number of parameters to two, namely 
X and E. While the hyperpolarizability is known to 
depend on the contributions from many excited states, it 
was proposed that when a quantum system is at the limit, 
only two excited states contribute. [l6[ Thus, Equation 
[2l1 applies, in principle, only at the limit. The second 
assumption is that X and E arc independent. Under this 
condition, f{E) and G{X) arc individually optimized and 
found to peak at unity for i? and X = sa 0.75. 

The solid black curve in Figure [S] represents 
ztf{E)G{X) for E fixed by the energy class. Also shown 
is ztf{E) as horizontal dashed curves. For each point on 
the plot, the number of dominant states is determined by 
the criteria that it contributes at least 25% to the total, 
and color coded red for three states, blue for 4 to 6 states, 
and green for seven or more states. First, we find that 
f{E), aside from a small number of outliers, determines 
the range of Pint values observed. As E gets smaller. 



more points in the distribution are closer to the limit. 

The magnitude of the hyperpolarizability is also found 
to correlate with the number of dominant states. For 
seven or more dominant state, the hyperpolarizability 
fills a narrow band centered on Pint = and with a 
width of about 0.1. The systems with four to six dom- 
inant states form a broader band while the three-level 
systems fill the range up to the values f{E)G{X). There 
are several interesting features of the data. First, there is 
no reason for the function f{E) to determine the upper 
bounds for systems with hypcrpolarizabilities far from 
the limit. However, it appears that the three-level re- 
sults accurately predicts the upper bounds. Thus, for a 
given energy class, the sum-rule constrained three-level 
model appears to yield the range of Pint for each value 
of X. 

All the points that fall on the curve f{E)G{X) for 
every energy class represents a the three-level model. In- 
terestingly, this is the case not only for systems with an 
energy class for Pint ~ 1, but even when the energy func- 
tion f{E) limits the hyperpolarizability to 1/10 of the 
fundamental limit. However, there are also three-state 
systems that are far below the f{E)G{X) curve. Thus, 
when a quantum system is described by three dominant 
states, this does not guarantee that the hyperpolarizabil- 
ity is at the fundamental limit even when the energy ratio 
and transition moments meet the criteria. This is most 
likely due to many other excited states that are each 
contributing less than 25% of the total, but in aggregate, 
decrease the value below the limit. 

It is interesting that a plot calculated from the sum- 
rulc-constraincd three-level model so accurately provides 
a bound for all of the observed data. In particular, in 
cases where the energy class yields a small value of f{E), 
why are there no points outside the figure-eight pattern? 
Furthermore, it is interesting that the hyperpolarizabili- 
ties all vanish when X = 0. While one may question if 
the three-level ansatz holds, our data here suggests that 
not only docs it hold as originally posed, but the three- 
level ansatz is more broadly applicable than just at the 
limit. Furthermore, the data suggests that indeed, X and 
E are independent, as previously postulated. 



C. Second Hyperpolarizability 

In this section we focus on the second hyperpolariz- 
ability. Figure [H] shows the distributions of jint for the 
energy class The range of the intrinsic second hy- 

perpolarizability spans the range —0.10 < "fint < 0.15. 
However, the three-level model clusters around the value 
of 7int = 0.05. Thus, the three-level model gives a result 
that is far less than systems with more levels. However, 
the four-level model peaks at the largest value of posi- 
tive jint- But, as we saw in the energy class e-', shown in 
Figure [H the three- level model peaks at the limit for the 
optimized energy spectrum. 

The three-level ansatz when applied to the off-resonant 
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FIG. 6. The distribution of the second hyperpolarizabihty for 
energy class j""^. 

second order hyperpolarizability, jmax yields [39| . 
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where 

G^{E, X)='i-5{E-lf{E + 1) {E^ + E+ l) X"^ 
-2 {E^ - 1) E^X^ - {E^ + E + 3) E^. 

(23) 

The maximum value of G^{E,X) is 4 which is obtained 
for the extreme case E ~ 0, leading to the fundamental 
limit of second hyperpolarizability, jmax- The intrinsic 
value of any arbitrary system is found by dividing Equa- 
tion [22] by jmax, defined in Equation [2l 
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In analogy to /3, we are interested in the dependence 
of G-y{E,X) on E. Therefore we introduce f{E) = 
G{E, XoiE))/4 where Xo{E) is the value of X that max- 
imizes G {E, Xo{E)) for fixed E. We find Xo{E) by solv- 
ing the equation dG/dX = for X{E), yielding 



X^q'\e)^0 and X^^\e)^±- ^^'^ 



(25) 



It is straightforward to show that the maximum value 
of G^ is obtained 
Equation [53] gives, 



of G-y is obtained for X^ {E). Introducing X^ '{E) in 
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FIG. 7. Distribution of -yint with energy for several differ- 
ent energy functions. The vertical dashed lines represent the 
energy function f'y{E). 



For example, for the energy class Ej = j ^ where E = 
E10/E20 = 0.843, we find f^{E) = 0.22. The vertical 
dashed lines in Fig. [7] represent f-y{E) values. 

Figure [7] summarizes the results of the second hyperpo- 
larizability calculations for the various energy classes. It 
is important to keep in mind that all the distributions for 
7mt are plotted on a log scale. As the energy spectrum 
becomes more spread out, the range of 7i„t increases, 
where the range of ^i„t of the three-level model is smaller 
than the range in "fint with 4 or more states. However, as 
the energy spectrum becomes more favorable, the three- 
state plot approaches the same range of jint values as 
the other plots with a peak at "fint approaching unity. 
In all simulations, the fundamental limit is not exceeded, 
suggesting that the three-level ansatz gives the correct 
limit. 

Table [ll] gives a sampling of the largest jint values, 
the dominant levels and their contribution to the intrin- 
sic second hyperpolarizability for the three energy classes 
Ej = exp(j), Ej = j and Ej ~ for a 20-state model. 
To find the contribution of pairs of states to 7i„t, we 
use the missing state analysis In this method, we 

calculate "fmt in the absence of the pairs of states i and 

The fractional difference be- 



(26) j, which is called 7, 



missing 
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TABLE II. Summary of the properties of the second hyperpo- 
larizabihty for a sampUng of energy functions when the second 
hyperpolarizability is small or near the limit. 



Function 


lint 


Levels Dominant States 




0.443 


5 1,2,3,4 


j 


0.92 


4 1,2,4 


exp(j) 


0.973 


3 1,2 



tween ^™*^'**"9 ^nd 7i„t describes the joint contribution 
of states i and j to ^int ■ The smaller the value of ry™^'^™s ^ 
the larger the contribution of states i and j to the second 
hyperpolarizability. [23] For example, for Ej — exp(j). 



^rmssirig ^j^^ order of 10 ^ and 



7]^ 4 IS the sec- 



ond most significant term, which is of the order of 10"^ 
- indicating that the pair (1, 2) contributes to the second 
hyperpolarizability 100 times more than the pair (1,4). 

For Ej = 5 states contribute significantly to 7 

with (2, 3) having the largest contribution. The energy 
class Ej = j has a larger value of E than Ej = j^^, but 
only four states dominate. For Ej = exp(j), for which 
the second hyperpolarizability is closest to the limit of 
the three cases, only three states (ground state and two 
excited states) contribute to 7, as predicted by three level 
ansatz. 

Figure|S]illustrates the distribution of 7i„t as a function 
of X for a 20-state model. While the Monte Carlo simu- 
lations generate all possible X values, the resultant jint 
range depends on the energy class. When E is large, e.g. 
Ej = and Ej = the largest intrinsic second hy- 
perpolarizabilities do not exceed 0.3 while for classes such 
as Ej = p ; they approach the limit. The dashed lines 
in Figure [8] are the maximum attainable values of 7i„t 
in the sum-rule restricted three-level model, i.e. f-y{E). 
The numerical values of f-y {E) are given in Figure [71 

When the system is not close to the fundamental limit, 
there are 7i„t values that exceed /-y [E) as represented by 
the vertical dashed lines in Figure [71 indicating that the 
sum-rule restricted three-level model does not necessar- 
ily produce the largest 7i„t values when the system is 
not optimized. Otherwise, when — > 0, the three-level 
model leads to 7,;„t — > 1 as is the ease for Ej = in 
Figure [3 

The curved solid lines in Figure [H plot {X, E) /4 
for the value of E that is defined for each energy class, 
indicating the largest 7i„t values that can be generated 
by each energy class when they are restricted to a sum- 
rule-constrained three-level model. When the three-level 
model is not optimized, i.e. G{X, E) 7^ 4, the distri- 
bution of 7i„t exceeds G.y{X, E)/4. However, as the en- 
ergy parameter E gets smaller, G.y{X, E)/A increases and 
fewer points in the distribution fall outside the energy 
function. When E = {), all points in the distribution fall 
below G^{X, 0)/4, which defines the fundamental limit of 
the hyperpolarizability. This behavior supports the va- 
lidity of three-level ansatz in calculating the fundamental 
limits of the second hyperpolarizability. 
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FIG. 8. 7i„t as a function of X for several energy classes 
using a 20-state model. Dashed lines represent F^,{E). The 
black solid curve specifies the limit defined by the sum-rule- 
constrained three-level model. 



D. The Gap 

The gap between the experimental results and the fun- 
damental limit for the firstjl^ and second]!^ order hy- 
perpolarizabilities have been extensively discussed in the 
literature in which Monte Carlo simulations suggest that 
the gap might be due to the unfavorable arrangement of 
the excited state energies. 

Figure [9l shows the distribution of the energy class 
and 1 which approximately represents the energy 
spacing of atoms and molecules. Most of the hyperpolar- 
izabilities fall in the range —0.1 < (imt < +0.1 for class 
j^^ and —0.1 < /3int < +0.1 for class j^^. Thus, since 
most organic molecules fall in this range of energy class, 
one would expect that the very best molecules would fall 
below lAntI « 0.2. 

Assuming that the transition moments vary more be- 
tween molecules than does the energy spacing. [6l[37| the 
statistically most likely observed hyperpolarizabilities are 
those near the peak in the distribution function. The 
energy classes and both have a double-peaked 
structure in the range —0.03 < /3i„t < +0.03 for the sum- 
rule constrained three-level model. When more states 
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FIG. 9. The distribution of /3i„t for the energy class j ^ (top) 
and j"'' (bottom). 



contribute, the double-peaked behavior is not observed, 
but the width of the peak remains at 0.03. Thus, based 
on statistics alone, one would expect most molecules to 
have an intrinsic hyperpolarizability of less than 0.03 - 
precisely the value of the observed gap. 

At the heart of the gap may be the nature of coulomb 
forces, which lead to an energy spacing that is not con- 
ducive to optimizing the nonlinear-optical response un- 
less the effective potential energy function is interrupted 
in a way to change the energy spacing of the lowest- 
energy eigenstates. 

In the case of the second hyperpolarizability, the 
spread in jint is larger than for Pint- Furthermore, the 
distribution is flatter so it is statically more likely to find 
molecules with larger intrinsic second hyperpolarizabil- 
ity. For the energy class for example, the distribu- 
tion narrows as the number of states in increased, but 
in all cases, the range is at least —0.07 < 7mt < 0.2. 
This is consistent with reports of jint ~ 0.2 by May and 
coworkers, [loj 

Finally, in numerical simulations in which the po- 
tential energy function is varied to optimize Pint, the 



largest values observed are pint ~ 0.708. [IJ, [ij, 122, l23| 
These simulations show that at the limit, potentials share 
certain universal properties, such as the energy ratio 
E10/E20 ~ 0.48. This corresponds to the energy class 
J, which shows a sharp drop-off beyond Pint = 0.707. 
Thus, it appears that it may not be possible to reach the 
fundamental limit by varying the potential alone in ID 
molecule. 



IV. CONCLUSION 

Classifying Monte Carlo simulations using an energy 
spectrum function resolves several long-standing ques- 
tions. First, our work shows the centrality of energy 
spacing in determining the intrinsic nonlinear response. 
While a broad range of transition moments are observed 
in atoms and molecules, the energy spacing - as charac- 
terized by the energy parameter, E, varies little between 
systems. Indeed, the importance of the energy parameter 
in attaining larger hyperpolarizabilities has been demon- 
strated in several experimental studies. [gI. Ist} 

The factor of 20-30 gap between the best molecules 
and the fundamental limit of the hyperpolarizability can 
be understood statistically. When the energy spacing is 
characterized by = 0, then a random sampling of quan- 
tum systems shows an approximately flat distribution of 
hyperpolarizabilities from Pint = to Pint = 1- However, 
when the energy function is i? « 0.8, as it is in typi- 
cal molecules, a random sampling of transition moments 
would yield hyperpolarizabilities that fall a factor of 20 
to 30 below the limit. The largest possible hyperpolariz- 
ability in this case is about 0.1 to 0.2. Unless quantum 
systems are identified that have a more favorable spac- 
ing, the intrinsic hyperpolarizability can be improved by 
no more than about a factor of 5. 

The second-order hyperpolarizability, on the other 
hand, shows a similar behavior but with less sensitiv- 
ity to the energy parameter. For E « 0.8, jint can be 
as large as about 0.2. This is consistent with the iden- 
tification of small molecules with intrinsic second hyper- 
polarizabilities in this range. [l^l Thus, when viewed 
in terms of energy classes, our Monte Carlo simulations 
explain the observed gap between the best molecules and 
the fundamental limit. 

Studies aimed at finding the optimum potential energy 
functions that maximize the intrinsic hyperpolarizabil- 
ity find many potential functions that represent a local 
maximum of /Jmt-O [iB] Quantum systems that are so 
optimize all share certain universal properties, such as 
E w 0.48 and Pint ~ 0.708. In our present studies, we 
find that the energy class Ej = j, with E = 0.5 is re- 
stricted to Pint < 0.708. As such, our work suggests 
that solutions to the Schrodinger Equation with a poten- 
tial energy function demands that the energy parameter 
E > 0.48. Thus, in the process of optimizing the poten- 
tial energy function, the energy parameter is optimized 
to its minimum possible value. This is consistent with the 
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suggestion that systems with more exotic Hamiltonians 
may be required to attain the fundamental limit. [2^ 

All of our simulations have shown the validity of the 
three-level ansatz - that is, at the fundamental limit, the 
system is represented by a sum-rule constrained three- 
level model with all of the oscillator strength shared by 
at most two states. We find no instances where more than 
three states contribute at the limit, though a three-state 
model may describe a system that is far from the limits. 
Also, since no values are found to exceed the limit, the 
theoretical foundations of the limit calculations appear 
to be on solid ground. 

More interestingly, the sum-rule constrained three- 
level model appears to be more broadly applicable than 
just at the limits. Even when the quantum system's hy- 
pcrpolarizability is far from the limit, the energy func- 
tion f{E) defines the limits of the nonlinear response, 
as one would expect if the functions f{E) and G{X) - 
which define the intrinsic hyperpolarizability according 



to (3int = f{E)G{X) - are independent. Similarly, G(X) 
also provides a constraints, such as forcing Pint to be 
small when X is small. In all cases, we find that as the 
number of states that contribute to the hyperpolarizabil- 
ity increases, the intrinsic hyperpolarizability is limited 
to narrower and narrower bands. Similar behavior has 
been observed for the second hyperpolarizability. 

Monte Carlo calculation using the energy classification 
scheme have bridged the divide between Monte Carlo 
simulations and potential energy optimization studies. 
The power of the Monte-Carlo technique lies in the fact 
that all possible Hilbert spaces are probed, leading to 
very broad and fundamental relationships. Using energy 
classifications allows the parameter space to be reduced 
to subsets that describe atoms and molecules. Future re- 
finements may lead to more specific design guidelines for 
making improved molecules for a variety of applications. 
The potential for discovering new fundamental science 
with this approach is of equal importance. 
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